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1. Technical Field 

The present invention relates to a method and algorithm for performing a sequence of 
integer divides without propagation of truncation error. 

2. Related Art 

Real-time video systems requiring large network bandwidth generally have their video 
signals compressed so that the video signals may be efficiently transmitted from source to 
destination. An example of an emerging video compression standard is the Moving Picture 
Experts Group (MPEG) standard. Under MPEG, video frames are initially encoded (i.e., 
compressed) for efficient transmission, placed in a buffer, and subsequently decoded (i.e., 
uncompressed) for viewing. 

During the encoding and decoding of video signals, a Video Buffering Verifier (VBV) 
buffer is dynamically filled with binary bits of encoded video data at a variable rate (in 
bits/frame), and the video data is subsequently removed from the buffer at a constant rate (in 
bits/frame) for decoding purposes. The constant bits/frame removal rate is referred to herein as 
an average bits/frame (BA), which is computed as 



Background of the Invention 



BA - BR/FR 



(1) 
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wherein BR is a bit rate in bits/second and FR is a frame rate in frames/sec. If integer arithmetic 
is used by buffer management software to compute BA, then Equation (1) should be recast into 
an integer format. For example if FR=29.97 frames/sec, then 
BA = BR/29.97 

= BR/(30-.03) 

= BR/(30-30/1000) 

= BR/(30(1-1/1'000)) (2) 
Expanding (1-1/1000) in a Taylor series, 

(1-1/1000)" 1 = 1 + 1/1000 + terms of second order and higher in 1/1000 (3) 
Thus to first order, substitution of Equation (3) into Equation (2) yields: 

BA = (BR + BR/1000)/30 (4) 
Equation (4) is a representation of Equation (1) to first order of 1/1000. Unfortunately, integer 
division by 1000 and by 30 in Equation (4) causes truncation error, which results in a smaller 
computed value of BA than is the "true" value of BA. The "true" value of BA is the constant 
number of bits/frame physically removed, while the smaller computed value of BA is the 
bits/frame that the buffer management software tracks as being removed based on Equation (4). 
For example, if BR=29970 bits/sec, then using Equation (4) with floating point arithmetic yields 
a "true" value of 1000 bits/frame (actually 999.999 bits/frame) for BA, but using Equation (4) 
with integer divides by 1000 and 30 yields a smaller computed of 999 bits/frame for BA. 
Accordingly with integer divides, the buffer management software would account for removal of 
video data at 999 bits/frame from the VBV buffer, while in actuality 1000 bits/frame is 
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physically removed from the VBV buffer. If B represents the number of bits stored in the VBV 
buffer at any given time, then B will be computed as a smaller value B c by the buffer 
management software than the true value B T . As video frames are processed, B T - B c will grow 
in magnitude and thus pose a risk that eventually B c will be so small in comparison with B T that 
the buffer management software will erroneously attempt to extract more bits from the VBV 
buffer than is actually in the VBV buffer. In accordance with the MPEG-2 standards, such an 
erroneous attempt would cause a VBV buffer violation labeled as "buffer overflow." 

Thus, there is a need to avoid buffer overflow during encoding and decoding of video 
frames of a video signal. 

Summary of the Invention 

The present invention provides a method for computing an average bits/frame (BA) for 
frames extracted from a buffer used for video encoding and decoding, each said frame having a 
same number of fields, said BA equal to (BR + BR1/J1)/J2, said BR1, Jl, and J2 each a positive 
integer, said BR a bit rate in bits/sec, said BR1/BR a positive integer, said method comprising: 

determining BR1, Jl, and 32 such that J2/(1+(BR1/BR)/J1) as evaluated in floating point 
is approximately equal to FR, said FR a frame rate in frames/sec; 

calculating a quotient Ql and remainder Rl from integer division of BR1 by Jl ; 

calculating a quotient Q2 and remainder R2 from integer division of (BR+Q1) by J2; 

initializing to zero accumulators Al and A2; and 

executing N iterations, wherein N > 1, and wherein executing each iteration includes: 
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adding Rl to Al; 

if Al ^ Jl , then adding 1 to A2 and decrementing Al by Jl ; 
setting BA=Q2 and adding R2 to A2; 
if A2^J2, then adding 1 to BA and decrementing A2 by J2. 
The present invention provides a computer code that computes an average bits/frame 
(BA) for frames extracted from a buffer used for video encoding and decoding, each said frame 
having a same number of fields, said BA equal to (BR + BR1/J1)/J2, said BR1, Jl, and J2 each a 
positive integer, said BR a bit rate in bits/sec, said BR1/BR a positive integer, said computer 
code including an algorithm programmed to: 
1(F: determine BR1, Jl, and J2 such that J2/(1+(BR1/BR)/J1) as evaluated in floating point is 

approximately equal to FR, said FR a frame rate in frames/sec; 

calculate a quotient Ql and remainder Rl from integer division of BR1 by Jl; 
calculate a quotient Q2 and remainder R2 from integer division of (BR+Q1) by J2; 
2 initialize to zero accumulators Al and A2; and 

La. 
[-— 

IS execute N iterations, wherein N > 1, and wherein to execute each iteration includes: 

to addRl to Al; 

if Al ^ Jl , then to add 1 to A2 and to decrement Al by Jl ; 
to set BA=Q2 and to add R2 to A2; and 
if A2^ J2, then to add 1 to BA and to decrement A2 by J2. 
20 The present invention provides a method of computing an average bits/frame (BA) for 

frames extracted from a buffer used for video encoding and decoding, each said frame having a 
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variable number of fields, comprising: 

defining BA1 as an average bits/frame for a two-field frame, said BA1 equal to (BR + 
BR1/J1)/J2, said BR1, Jl, and 32 each a positive integer, said BR a bit rate in bits/sec, said 
BR1/BR a positive integer; 
5 defining BA2 as an average bits/frame for a one-field frame, said BA2 equal to (BR + 

BR1/J1)/(2*J2); 

determining BR1, Jl, and J2 such that J2/(l-f(BRl/BR)/Jl) as evaluated in floating point 
is approximately equal to FR, said FR a frame rate in frames/sec; 

O 

y3 calculating a quotient Ql and remainder Rl from integer division BR1/J1 ; 

l&n calculating a quotient Q2 and remainder R2 from integer division (BR+Q1)/J2; 
UJ 

In 

~ calculating a quotient Q3 and remainder R3 from integer division (BR+Q1)/(2*J2); 

CO 

3 initializing to zero accumulators A 1 , A2, B 1 , and B2; 

ill executing N iterations, wherein N > 1, said executing iteration n of N relating to 

L-JL 

^ extracting a frame n from the buffer, said executing of iteration n including: 
1 5 calculating B A 1 , including : 

adding Rl to Al; 

if Al £ Jl then adding 1 to A2 and decrementing Al by Jl ; 
setting BA1=Q2 and adding R2 to A2; 
if A2^J2, then adding 1 to BA1 and decrementing A2 by J2; 
20 determining a number of fields F n comprised by the frame n; 

if F n is even then setting BA2=0 else calculating BA2 including: 
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adding Rl toBl; 

if B 1 ^ Jl ? then adding 1 to B2 and decrementing Bl by Jl; 
setting BA2=Q3 and adding R3 to B2; 

if B2^(2*J2), then adding 1 to BA2 and decrementing B2 by (2*J2); 
computing BA=(F n /2)*BAl + BA2, said (F n /2) computed by integer division. 
The present invention provides a computer code that computes an average bits/frame 
(BA) for frames extracted from a buffer used for video encoding and decoding, each said frame 
having a variable number of fields, said BA a function of BA1 and BA2, said BA1 defined as an 
average bits/frame for a two-field frame, said BA1 equal to (BR + BR1/J1)/J2, said BR1, Jl 5 and 



1 (Jfl J2 each a positive integer, said BR a bit rate in bits/sec, said BR1/BR a positive integer, said 

UJ 

LH 

q B A2 defined as an average bits/frame for a one-field frame, said BA2 equal to (BR + 

s BR1/J1)/(2*J2), said computer code including an algorithm, said algorithm programmed to: 
O 

Ul determine BR1, Jl, and J2 such that J2/(R(BR1/BR)/J1) as evaluated in floating point is 

^ approximately equal to FR, said FR a frame rate in frames/sec; 

O 

Us 

1 § calculate a quotient Ql and remainder Rl from integer division BR1/J1 ; 

calculate a quotient Q2 and remainder R2 from integer division (BR+Q1)/J2; 
calculate a quotient Q3 and remainder R3 from integer division (BR+Q1)/(2*J2); 
initialize to zero accumulators Al, A2, Bl, and B2; 

execute N iterations, wherein N > 1, said iteration n of N relating to extracting a frame n 
20 from the buffer, wherein to execute iteration n includes: 

to calculate BA1, including: 
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to addRl to Al; 

if Al £ Jl then to add 1 to A2 and to decrement Al by Jl ; 
to set BA1=Q2 and to add R2 to A2; 
if A2^J2, then to add 1 to BA1 and to decrement A2 by J2; 
to determine a number of fields F n comprised by the frame n; 
if F n is even then to set BA2=0 else to calculate BA2 including: 
to addRl toBl; 

if Bl^Jl, then to add 1 to B2 and to decrement Bl by Jl; 
to set BA2=Q3 and to add R3 to B2; 
l(|y to compute BA=(F n /2)*BAl + BA2, said (F n /2) computed by integer division. 

m 

Q The present invention provides a method for computing Z, said Z = E n Z n? said E n 

m 

* denoting a summation over n from 1 to N, said N a positive integer of at least 1, said Z n = X n /Y, 

m 

M said X n = (I ln /J,)M ln + (I 2n /J 2 )M 2n + ... + 

M 1 

° (WJi<;)M Kn , said Y and said I kn , J k , M kn (k=l, 2, ... 5 K) each a positive integer, said K a positive 
1 5 integer of at least 1 ? said method comprising: 

setting Z=0 5 B=0, and A k =0 for k=l, 2 3 K; 

executing N iterations, said executing of iteration n of N including: 

calculating a quotient Q kn and a remainder R kn from integer division I kn /J k for k=l 5 
2, K; 

20 calculating X n = E k [Q kn M kn ] as summed over k from 1 to K; 

adding R kn M kn to A k for k=I, 2, K; 
END920010005US1 7 



for k = 1 , 2, K, if A k ^ J k5 then adding 1 to B and decrementing A k by J k ; 
if Y* 1 then calculating a quotient Q n and a remainder R n from integer division 
X n /Y, else setting Q n = X n and R n = 0; 

setting Z n = Q n and adding R n to B; 

if B ^Y, then calculating Z n = Z n + 1 and decrementing B by Y; 
adding Z n to Z. 

The present invention provides a computer code that computes Z, said Z = E n Z n? said 
E n denoting a summation over n from 1 to N, said N a positive integer of at least 1 , said Z n = 
X n /Y, said X n = (I In /J,)M ]n + (I 2n /J 2 )M 2n + ... + (I Kn /J K )M Kn , said Y and said I kn , J k , M kn (k=l, 2, .... 
K) each a positive integer, said K a positive integer of at least 1 , said computer code including an 
algorithm, said algorithm programmed to: 

set Z=0, B=0, and A k =0 for k=l, 2, K; 

execute N iterations, wherein to execute iteration n of N includes: 

to calculate a quotient Q kn and a remainder R kn from integer division I kn /J k for k=l, 
2, K; 

to calculate X n = E k [Q kn M kn ] as summed over k from 1 to K; 
to add R kn M kn to A k for k=l, 2, K; 

for k = 1, 2, K, if A k ^ J k , then to add 1 to B and to decrement A k by J k ; 
if Y* 1 then to calculate a quotient Q n and a remainder R n from integer division 
X/y, else to set Q n = X n and R n = 0; 

to set Z n = Q n and to add R n to B; 
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if B ^ Y, then to calculate Z n = Z n + 1 and to decrement B by Y; 
to add Z n to Z. 

The present invention avoids buffer overflow during encoding and decoding of video 
frames of a video signal. 

Brief Description of the Drawings 

FIG. 1 is a flow chart showing an iterative calculation of extracted bits/frame 
(BA_FINAL) of video processing in terms of a single-picture bits/frame (BA_FINAL1) and half- 
picture bits/frame single-picture bits/frame (BA_FINAL2), in accordance with embodiments of 
the present invention. 

FIG. 2 is a flow chart for calculation of the single-picture bits/frame (BA_FINAL1) of 
FIG. 1, in accordance with embodiments of the present invention. 

FIG. 3 is a flow chart for calculation of the half-picture bits/frame (BA_FINAL2) of 
FIG. 1, in accordance with embodiments of the present invention. 

FIG. 4 is a table illustrating iterations of FIG. 1 for the BA_FINAL1 calculation of FIG. 2 
for an example, in accordance with embodiments of the present invention. 

FIG. 5 is a table illustrating iterations of FIG. 1 for the BA_FINAL2 calculation of FIG. 3 
for an example, in accordance with embodiments of the present invention. 

FIG. 6 is a flow chart for an iterative calculation with integer division, in accordance with 
embodiments of the present invention. 

FIG. 7 is a flow chart for an iterative summation calculation with integer division, in 
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accordance with embodiments of the present invention. 

FIG. 8 is a table illustrating iterations of FIG. 7 for an example, in accordance with 
embodiments of the present invention. 

FIG. 9 illustrates a computer system, in accordance with embodiments of the present 
invention. 

Detailed Description of the Invention 

Video compression under MPEG-2 includes three types of compression: spatial 
compression, temporal compression, and 3/2 compression. Spatial compression compresses data 
within a picture such as by use of discrete cosine transformation (DCT) coefficients to account 
for spatial redundancy within the picture. Temporal compression compresses data between 
successive pictures such as through motion compensation (e.g., describing a picture in terms of 
vectors that relate portions of the picture to corresponding portions of the previous picture). 3/2 
compression compresses V/ 2 pictures into 1 picture as will be explained as follows. 

Under MPEG-2, there is a distinction between pictures, fields, and frames. A picture 
comprises a physical screen of binary bits. Under MPEG-2, a picture consists of two interlacing 
fields, namely an upper field that includes bits on an upper portion of the picture and a lower 
field that includes bits on a lower portion of the picture. Each such field may be thought of as a 
half-picture. A frame is collection of fields that is stored in a Video Buffering Verifier (VBV) 
buffer as a unit. Under MPEG-2, a frame has either 2 fields or 3 fields and is called a 2-field 
frame or a 3-field frame, respectively. An example of a 2-field frame is the upper and lower 
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fields of a single picture. An example of a 3-field frame is the upper and lower fields of a first 
picture followed by the upper field of a second (i.e., next) picture. With a 3-field frame having 
fields sequentially denoted as fields 1, 2, and 3, the field 3 may be indistinguishable to a human 
eye from the field 1 . Thus with the previous example of a 3-field frame, the upper field of the 
second picture is indistinguishable to the human eye from the upper field of the first picture. 
Accordingly, field 3 of the 3-field frame is represented by a "repeat flag" that denotes repetition 
of the bits of field 1. Since the 3-field frame includes only 2 fields of bits and the repeat flag, the 
3-field frame represents a compression known as the "3/2 compression." 

During video processing of the VBV buffer for subsequent decoding, BA bits/frame are 
removed from the VBV buffer. Each such removal of BA bits from the VBV buffer constitutes 
an "iteration" in a sequence of such removals. Upon such removal in an iteration, the buffer 
management software must compute B A corresponding to the average number of bits/frame 
removed and keep track of the number of bits present in the VBV buffer at the end of each 
iteration. Equation (4) may be viewed as a calculation for BA corresponding to a 2-field frame 
of a single picture. Thus if a 2-field frame is processed during an iteration, then Equation (4) 
may be used for calculating BA. But if a 3 -field frame is processed during an iteration, then BA 
corresponding to 3 fields, or VA pictures, must be calculated. Thus a 3-field frame value of BA 
is equal to BA of a single picture + half-BA of a single picture, which represents BA for 1 Vi 
pictures. Accordingly, calculation of BA for a given iteration depends on whether the frame of 
extracted bits for the iteration is a 2-field frame or a 3-field frame, as depicted in FIGS. 1-3, in 
accordance with embodiments of the present invention. FIGS. 1-3 are based on Equation (4) 
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with modifications in accordance with the present invention. 

FIG. 1 is a flow chart illustrating looping through iterations. In each iteration, the total 

bits/frame (BA_FINAL) is calculated to represent the average number of bits/frame that is 

extracted from the buffer in each iteration. As seen in block 20, BA_FINAL is a sum of terms 

5 BAJFINALl and BA_FINAL2. The term BA_FINAL1 represents BA for a 2-fields (i.e., full 

picture) and is calculated as shown in FIG. 2. The term BA_FINAL2 represents BA for a 1 -field 

frame (i.e., half-picture) and is calculated as shown in FIG. 3. For a 2-field frame, 

BA_FINAL2=0 and BA_FINAL = BA_FINAL1. For a 3-field frame, BA FINAL is a sum of 

|! BA for one picture (i.e., BA_FINAL1) and BA for a half-picture (i.e., BA_FINAL2), which in 
in 

lGbj composite represents 1 X A pictures, 
g Block 10 of FIG. 1 includes initializations for the BA_FINAL1 calculations of FIG. 2 and 

q initializations for the BAFFIN AL2 calculations of FIG. 3. The initialization block 10 will be 

m 

H described infra in conjunction with FIGS. 2 and 3. Blocks 12, 14/16 or 14/18, and 20 are 

K 

P executed within a single iteration. Thus, the blocks 12, 14/16 or 14/18, and 20, together with a 
15 return path 22, defines an iteration loop. Block 12 calculates BA_FINAL1 as described in FIG. 

2. Decision block 14 asks whether the frame being processed for extraction is a 3 -field frame. If 
NO, then BA_FINAL2 is set equal to 0 (for a 2-field frame) as shown in block 16. If YES, then 
block 18 is executed, and block 18 calculates BA_FINAL2 (for a 3-field frame) as described in 
FIG. 3. Block 20 sums BA_FINAL1 and BA FINAL2 to calculate BA_FINAL for either a 2- 
20 field frame or a 3 -field frame. The path 22 effectuates a transition between successive iterations. 

The flow chart of FIG. 1 may be modified in any manner that is logically equivalent to 
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FIG. 1 as shown herein, as would be understood by one of ordinary skill in the art. For example, 
blocks 16 and 18 could be replaced by the following logic. If NO is the answer to the query in 
decision block 14, then BA_FINAL =BA_FINAL1 is executed. If YES is the answer to the 
query in decision block 14, then BA_FINAL2 is calculated followed by execution of BA FINAL 
= BA_FINAL1 + BA_FINAL2. 

FIG. 2 is a flow chart for calculation of BA_FINAL1 of the block 12 of FIG. 1. Thus the 
flow chart of FIG. 2 is within the iterations loop shown in FIG. 1. In FIG. 2, BA_FINAL1 is 
calculated, as shown in FIG. 2, as a function of BA1 : 



wherein BR is bit rate in bits/second expressed as a positive integer, wherein FR_A is a frame 
rate in frames/sec (e.g., frames/sec), and wherein the divisions by 1000 and FR_A are integer 
divisions with truncation of fractional remainders. It is assumed herein that "integer division" 
generally results in truncation of fractional remainders. 

The initialization 10 of FIG. 1 comprises the following initializations for the 
BA FINAL 1 calculation of FIG. 2: 



BA1 = (BR+BR/1000)/FR_A 



(5) 



adj_t = BR/1000 



(6) 



adt_rem = BR - (1000*adj_t) 



(7) 



ba_t = (BR+adj_t)/FR_A 



(8) 



ba_rem = (BR+adj_t) - (FR_A*ba_t) 



(9) 



adj_rem_accum = 0 



(10) 



ba rem accum = 0 



in) 
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wherein the divisions by 1000 and FR_A in Equations (6) and (8), respectively, are by integer 
division with truncation. For example, if BR=26320, BA1 would equal 878.21 in a floating 
point implementation of Equation (5). With integer division, however, the initializations of 
Equations (6)-(9) yield: 
5 adj_t=26 

adj_rem=320 
ba_t=878 
_ ba__rem-6 

^ Note that baj in Equation (8) is the truncated value of BA1 of Equation (5), as illustrated in the 

m 

IOjJ preceding example such that the ba_t=878 and BA1 (floating point)= 878.2 1 . As shown in the 

5 block 66 of FIG. 2, BA_FINAL1 is equal to ba_t unless a parameter e adder is equal to 1 . The 

CO 

q parameter e_adder offsets the aforementioned truncation as will now be demonstrated. 

in 

The remainders adj_rem and ba_rem are lost through truncation resulting from integer 
U division in Equations (6) and (8), respectively. To compensate for this truncation, the algorithm 

15 of FIG. 2 saves and accumulates remainders adj_rem and ba_rem in accumulators 

adj_rem_accum and ba_rem_accum, respectively, in each iteration of FIG. 1, as shown in blocks 
32 and 48, respectively, of FIG. 2. An "accumulator" is defined herein as a storage location for 
storing a cumulative quantity (i.e., a summation). If in an iteration of FIG. 1, ba_jem_accum 
accumulates to FR_A (e.g. 30) or greater, then the accumulator ba_rem_accum is decremented 

20 by FR_A and BA_FINAL1 is incremented by 1 (i.e., e_adder=l in FIG. 2), as shown in blocks 
52, 58, and 62. If ba_rem_accum is less than FR_A, then e_adder is set to zero (see blocks 52 
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and 56) and BA^FINALl is not incremented by 1. Similarly, if in an iteration of FIG. 1, 
adj_rem_accum in FIG. 2 accumulates to 1000 or greater, then the accumulator adj_rem_accum 
is decremented by 1000 and the accumulator ba_rem_accum is incremented by 1 (i.e., 
ba_rem_adder=l in FIG. 2), as shown in blocks 34, 42, and 46. If adj_rem_accum is less than 
5 1000 in an iteration, then ba_rem_adder is set to zero (see blocks 34 and 38) and ba_rem_accum 
is not incremented by 1 . It should be noted from block 48 that ba__rem_adder is an increment to 
the accumulator ba_rem_accum. Thus the accumulators adj_rem_accum and ba_rem_accum 
compensate for the truncation losses in computing BA1 of Equation (5) by integer division of 

03 1 000 and FR_A, respectively. 

Ill 

1(W The flow chart of FIG. 2 may be modified in any manner that is logically equivalent to 

m 

gjl FIG. 2 as shown herein, as would be understood by one of ordinary skill in the art. For example, 

3 

O some or all of the intitializations of Equations (6)-(9), or mathematical equivalents thereof, could 
W 

M 1 alternatively be performed within the iteration loop of FIG. 2. Since adj _rem_accum and 

jp; ba_rem_accum are variables rather than constants, Equations (10) and (1 1) (or mathematical 
15 equivalents thereof) could each alternatively be performed within the iteration loop of FIG. 2 but 
only at the beginning of the first iteration. 

FIG. 4 illustrates for the previous example (i.e., BR=26320) values at the end of each of 
the first 1 1 iterations of: adj_rem_accum, ba_rem_adder, ba_rem_accum, e_adder, and 
BA FINALl. In FIG. 4, BA_FINAL1 has the truncated value of 878 at the end of iterations 1-4, 
20 6-9, and 11, but has a value of 879 in iterations 5 and 10. Hence, the algorithm of FIG. 2 

mathematically extracts an extra bit in iterations 5 and 10 to compensate for the truncations in 
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iterations 1-4 and 6-9. Thus, the algorithm of FIG. 2 prevents propagation of the truncation error 
by periodically simulating the effect of near-perfect integer division. 

FIG. 3 is a flow chart for calculation of BA_FINAL2 of the block 18 of FIG. 1 . Thus the 
flow chart of FIG. 3 is within the iterations loop shown in FIG. 1. In FIG. 3, BA_FINAL2 is 
5 calculated, as shown in FIG. 3, as a function of BA2: 

BA2 = (BR+BR/1000)/(2*FRA) (12) 
wherein the divisions by 1000 and 2*FR_A are integer divisions with truncation of fractional 
q remainders. 

CO The initialization 10 of FIG. 1 comprises, in addition to initializations of Equations (6)- 

1 (jfjl ( 1 1 ), the following initializations for the BA_FINAL2 calculation of FIG. 3 : 
Q 

m h_adj_t = BR/1000 (13) 

h_adj_rem = BR - (1000*h_adj_t) (14) 
h_ba_t - (BR+h_adj_t)/(2*FR_A) (1 5) 

h_ba_rem - (BR+h_adjJ) - (2*FR_A * h_ba_t) (16) 
1 5 h_adj__rem__accum = 0 (17) 

h_ba_rem_accum = 0 (18) 
wherein the divisions by 1000 and 2*FR_A in Equations (13) and (15), respectively, are by 
integer division with truncation. Note that h_adj_t and h_adj_rem of FIG. 3 are equal to adj__t 
and adj_rem of FIG. 2. For the preceding example with BR=26320, BA2 would equal 439.105 
20 in a floating point implementation of Equation (12). With integer division, however, the 
initializations of Equations (13)-(14) yield: 
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h_adj_t=26 
h_adj_rem=320 
h_ba_t=439 
h_ba_rem=6 

Note that h_ba_t in Equation (15) is the truncated value of BA2 of Equation (12), as illustrated in 
the preceding example such that the h_ba_t=439 and BA2 (floating point)= 439.105. 
Accordingly, B A_FINAL2 is equal to h_ba_t unless a parameter h_adder is equal to 1 , as shown 
in the block 166 of FIG. 3. The parameter h_adder offsets the aforementioned truncation as will 
now be demonstrated. 

The remainders h_adj_rem and h_ba_rem are lost through truncation resulting from 
integer division in Equations (13) and (15), respectively. To compensate for this truncation, the 
algorithm of FIG. 3 saves and accumulates remainders h_adj_rem and h_ba_rem in accumulators 
h_adj_rem_accum and h_ba_rem_accum, respectively, in each "YES" iteration of FIG. 1 (i.e., an 
iteration in which the answer to the question in the decision block 14 of FIG. 1 is "YES"), as 
shown in blocks 132 and 148, respectively, of FIG. 3. If in such a "YES" iteration of FIG. 1, 
h_ba_rem_accum accumulates to 2*FR_A (i.e. 60) or greater, then the accumulator 
h_ba_rem_accum is decremented by 2*FR_A and BAJFINAL2 is incremented by 1 (i.e., 
h_adder=l in FIG. 3), as shown in blocks 152, 158, and 162. If h_ba_rem_accum is less than 
2*FR_A, then h_adder is set to zero (see blocks 152 and 156) and BA_FINAL2 is not 
incremented by 1 . Similarly, if in such a "YES" iteration of FIG. 1, h_adj_rem_accum in FIG. 3 
accumulates to 1000 or greater, then the accumulator h_adj_rem__accum is decremented by 1000 
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and the accumulator h_ba_rem_accum is incremented by 1 (i.e., h_ba_rem_adder=l in FIG. 3), 
as shown in blocks 134, 142, and 146. If h_adj_rem_accum is less than 1000 in an iteration, then 
h_ba_rem_adder is set to zero (see blocks 134 and 138) and h_ba_rem_accum is not incremented 
by 1 . It should be noted from block 148 that h_ba_rem_adder is an increment to the accumulator 
h_ba_rem_accum. Thus the accumulators h_adj_rem_accum and h_ba_rem_accum compensate 
for the truncation losses in computing BA2 of Equation (12) by integer division of 1000 and 
2*FR_A, respectively. 

The flow chart of FIG. 3 may be modified in any manner that is logically equivalent to 

83 FIG. 3 as shown herein, as would be understood by one of ordinary skill in the art. For example, 

111 

i(W some or all of the intitializations of Equations (13)-(16), or mathematical equivalents thereof, 
could alternatively be performed within the iteration loop of FIG. 3. Since h_adj_rem_accum 
and h_ba_rem_accum are variables rather than constants, Equations (17) and (18) (or 
mathematical equivalents thereof) could each alternatively be performed within the iteration loop 
of FIG. 3 but only at the beginning of the first iteration. 
15 FIG. 5 illustrates for the previous example (i.e., BR=26320) values at the end of each of 

"YES" iterations (i.e., iterations 3, 4, 6, 9-11, 15, 18-19, and 22-23) of: h_adj_rem_accum, 
h_ba_rem_adder, h_ba_rem__accum, h_adder, and BA_FINAL2. In FIG. 5, BA_FINAL2 has the 
truncated value of 439 at the end of iterations 3, 4, 6, 9-11, 15, 18-19, and 23, but has a value of 
440 in iteration 22. Hence, the algorithm of FIG. 3 mathematically extracts an extra bit in 
20 iteration 22 to compensate for the truncations in iterations 3, 4, 6, 9-1 1, 15, and 18-19. Thus, the 
algorithm of FIG. 3 prevents propagation of the truncation error by periodically simulating the 
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effect of near-perfect integer division. 

While FIG. 1 implements 3/2 compression based on video encoding of either 2 
fields/frame or 3 fields/frame, the scope of the present invention includes a more general F/2 
compression, wherein F is the number of fields per frame such that F ^3. With F/2 compression, 
5 the question in decision box 14 of FIG. 1 is replaced by the following question: "Does this frame 
have F fields, wherein F ;>3?" Block 18 calculates BA_FINAL2 only if F is odd. Block 20 is 
changed to: 

BA_FINAL = (F/2) BA_FINAL1 + BA_FINAL2 (if F is odd) OR 
|j B AFFINAL = (F/2) B A_FINAL 1 (if F is even) 

1 <Xy wherein (F/2) is performed by integer division with truncation. For example, if F=3, which 



yi 

Q 

w 

3 



defines 3/2 compression, F/2=l in integer division so that BA_FINAL = BA_FINAL1 + 
BA_FINAL2 as shown in block 20 of FIG. 1. As another example, if F=4, which defines 4/2 



ji compression, F/2=2 in integer division so that BA_FINAL = 2*BA_FINAL1. While MPEG-2 
P currently limits F to 3, the present invention envisions inevitable advances in technology having 
15 improved compression such that F > 3. 

Equation (4) can be cast into a first different form by considering an example of 
FR=29.98 frames/sec in Equation (1). Using the same methodology that derived Equation (4), 
the following equation is derived: 

BA = (BR + (2*BR)/3000)/30 (19) 
20 Equation (4) can be cast into a second different form by considering an example of 

FR=29.94 frames/sec in Equation (1). Using the same methodology that derived Equation (4), 
END920010005US1 19 



the following equation is derived: 

BA = (BR + BR/500)/30 (20) 
Equation (4) can be cast into a third different form by considering an example of 
FR=19.97 frames/sec in Equation (1). Using the same methodology that derived Equation (4), 
the following equation is derived: 

B A = (BR + (3 *BR)/2000)/20 (2 1 ) 

Equations (4), (12), and (19)-(20) may be considered special cases of the following more 
general form: 

BA = (BR + BR1/J1)/J2 (22) 

m 

lCW wherein positive integers BR1, Jl, and 32 are determined such that J2/(1+(BR1/BR)/J1) as 

s 

J- evaluated in floating point is approximately equal to FR, and wherein BR1/BR is a positive 

p integer. BR1, Jl, and J2 may be so determined by, inter alia: being calculated, being received as 

m 

M* input, or by being hard-coded within an algorithm that implements Equation (22). Equations (4), 
Q 

r: (12), and (19)-(21) have the following common features: BR1/BR is a positive integer, Jl is a 
15 multiple of 10, and Jl > J2. BA as computed from Equation (22) is "approximately equal" to 

BR/FR (see Equation (1) for FR), which means that BA as computed from Equation (22) equals 
BR/FR to an extent that terms of second or higher order in a (i.e., the terms a 2 , a 3 ,...) are 
negligible in comparison with linear terms in a in the Taylor series for (1-cc)" 1 wherein a = 
(BR1/BR)/J1, as illustrated for ce=l/1000 in the derivation of Equation (4) from Equation (1). 
20 Note that ct=2/3000, 1/500, and 3/2000 in Equations (19), (20), and (21), respectively. Similarly, 
J2/(1+(BR1/BR)/J1) is "approximately equal" to FR. For a given FR and BR, there may be more 
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than one combination of BR1, Jl, and J2 such that J2/(1+(BR1/BR)/J1) is "approximately equal" 
to FR. 

Equation (22) is implemented as shown in the flow chart of FIG. 6, in accordance with 
embodiments of the present invention. The following initializations occur in block 70 of FIG. 6: 

a quotient Ql and a remainder Rl are calculated from integer division of BR1 by Jl; 

a quotient Q2 and a remainder R2 are calculated from integer division of (BR+Q1) by J2; 

and a first accumulator Al and a second accumulator A2 are each initialized to zero. 
Following the aforementioned initializations, N iterations are executed, wherein N > 1 . Each 
iteration includes blocks 72, 74, 76, 78, 80, 82, and 84, and return path 86. Block 84 executes a 
test that determines whether any more iterations remain to be executed. Executing each iteration 
includes: 

adding Rl to Al (block 72); 

if Al ^ Jl then adding 1 to A2 and decrementing Al by Jl ( blocks 74 and 76); 
calculating BA as equal to Q2 (block 78); 
adding R2 to A2 (block 78); and 

if A2^J2 then adding 1 to BA and and decrementing A2 by J2 ( blocks 80 and 82). 

Since Equation (22) encompasses Equations (5) and (12) as special cases, Equation (22) 
and the flow chart and algorithm of FIG. 6 could be used to implement the calculations of 
BA_FINAL1 and BA_FINAL2 of FIG. 2 and FIG. 3, respectively. Equivalently, Equation (22) 
and the flow chart and algorithm of FIG. 6 could be used to implement calculations of bits/frame 
(BA1) for two fields per frame (e.g., Equation (5)) and calculations of bits/frame (BA2) for one 
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field per frame (e.g., Equation (6)). For calculating BA1=BA in Equation (22), the positive 
integers BR1, Jl, and J2 are determined such that J2/(1+(BR1/BR)/J1) as evaluated in floating 
point is approximately equal to FR. For calculating BA2=BA in Equation (22) in a manner that 
is consistent with calculating BA1, Jl is replaced by 2* Jl . 

The flow chart of FIG. 6 may be modified in any manner that is logically equivalent to 
FIG. 6 as shown herein, as would be understood by one of ordinary skill in the art. For example, 
some or all of the intitializations block 70, or mathematical equivalents thereof, could 
alternatively be performed within the iteration loop of FIG. 6. Since accumulators Al and A2 
are variables rather than constants, initializations of A 1=0 and A2 =0 (or mathematical 
equivalents thereof) could each alternatively be performed within the iteration loop of FIG. 6 but 
only at the beginning of the first iteration. 

Equations (5) and (12), as well as Equation (22), which are computed with integer 
division by the algorithms described by FIGS. 1-3 and 6, may be expressed in a more general 
form as follows, in accordance with the present invention: 

Z n = X n /Y (23) 

X n = (I ln /J,)M In + (I 2n /J 2 )M 2n + ... + (I K n/JK)M K n (24) 
where "n" is an iteration index defining iterations analogous to the iterations described supra for 
FIG. 1 and FIG. 6. I kn , J k5 and M kn (k=l, 2, K; K ^ 1) are positive integers, and I kn /J k is 
performed by integer division with truncation. X n and Y are positive integers and X n /Y is 
performed by integer division with truncation. Additionally, Y and J k (k=l , 2, K) cannot be 
zero. As indicated in Equations (23)-(24), I kn and M kn (k=l, 2, K) are permitted to vary with 
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iteration index n. Y and J k (k=l, 2, K) are assumed to be constant and thus do not vary with 
iteration index n. Z n varies with iteration index n even if I kn and I kn do not vary with iteration 
index n, because the use of accumulators in accordance with the present invention causes Z n to 
increase at selected iterations in order to prevent propagation of error. 

The present invention could also be used to calculate the following summation: 
Z= £ n Z n (25) 
wherein the summation over n in Equation (25) is from 1 to N, and wherein N is the total number 
of iterations. 

Equations (23)-(24) reduce to Equation (5) if Y = FR_A, K = 2, I ln =BR, J,=l, M ln =l, 
I 2n =BR ? J 2 =1000, and M 2 „=l. Note that 3/2 compression, or more generally F/2 compression 
described supra, could be modeled in Equations (23)-(24) by having M kn =0 for those iterations in 
which the frame being processed has 2 fields, and M kn =l for those iterations in which the frame 
being processed has 3 fields, and suitable choices for Y, I k , J k , and M k . For example, 3/2 
compression, in combination with 2/2 compression, could be placed in the form of 

BA n = [BR + BR/1000 + (BR/2)*J n + (BR/2000)* J n ]/FR_A 
or BA n = [2*BR + BR/500 + BR*J n + (BR/1000)*J n ]/(2*FR_A) 
such that J n = 0 for iterations n having 2 fields/frame, and J n = 1 for iterations n having 3 
fields/frame. Note that the preceding equations for BA n conform to the form of Equations (23)- 
(24). 

Relative to FIGS. 1-3 and 6, and accompanying equations and text pertaining thereto as 
described supra, BR and FR may each be predetermined or received as input. 
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The present invention computes Z n in Equations (23)-(24) as described by FIG. 7. An 
accumulator B accumulates remainders resulting from integer division by Y in Equation (23). 
Accumulators A b A 25 ... 5 A K accumulate remainders resulting from integer division by J,, J 2s 
J K , respectively, in Equation (24). Block 200 in FIG. 7 initializes to zero the accumulators B and 
5 A k (k=l, 2, K) as well as the summation Z defined in Equation (25), as shown in Equations 

(23)-(24). Looping begins at block 210, and each loop iteration n includes blocks 210, 220, 230, 
240, 250, 260, 270, 280, 290, 300, and return path 310. As described by Equation (29), block 
210 calculates I kn /J k by integer division to yield an integer quotient Q kn and an integer remainder 

O 

g R kn such that R kn < J k for k=l , 2, K. In block 220 as described by Equation (30), X k is 

l(jfj computed by the summation E k5 from k=l to k=K, of Q kn M kn . Also in block 220, Equation (3 1 ) 

b I 

Q updates the accumulators A k for k=l, 2, K. Decision block 230 and block 240 are executed in 

3 sequence K times (i.e., blocks 230 and 240 for k=l, blocks 230 and 240 for k=2, blocks 230 
b 

'fj and 240 for k=K). Block 250 is not executed until blocks 230 and 240 are executed in sequence 
Q K times. Decision block 230 asks whether A k £ J k . If YES then B is incremented by 1 and A k is 
15 decremented by J k , as described by Equations (32) and (33), respectively, in block 240. If NO, 
block 240 is bypassed. 

As shown in Equation (34), block 250 calculates X n /Y by integer division to yield an 
integer quotient Q n and an integer remainder R n such that R n < Y. In block 260 as described by 
Equation (35), Z n is set equal to Q n . Also in block 260, Equation (36) updates the accumulator B. 
20 Decision block 270 asks whether B ^ Y. If YES then Z n is incremented by 1 and B is 

decremented by Y, as described by Equations (37) and (38), respectively, in block 280. If NO, 
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block 280 is bypassed. In block 290, Equation (39) updates the calculation of Z in accordance 
with Equation (25) supra. 

Decision block 300 asks whether there are more iterations to execute. If YES, the return 
path 310 directs the processing back to block 210 to begin the next loop iteration. If NO, the 
processing breaks the loop and executes post-processing in block 320. The post-processing 
block may adjust the calculation of Z to reflect the fact that after completion of all iterations, the 
accumulators B and A k (k=l, 2, K) may have non-zero integer contents. Thus Z may be 
adjusted in the post processing block 320 according to: 

Z = Z4-[B + E k (A k /J k )]/Y (40) 
wherein the summation over k in Equation (40) is from 1 to K. All calculations in Equation (40) 
should be performed in floating point. Thus, after Equation (40) is implemented, Z will be a 
floating point, or decimal, number. 

The flow chart of FIG. 7 may be modified in any manner that is logically equivalent to 
FIG. 7 as shown herein, as would be understood by one of ordinary skill in the art. For example, 
if any J k > = 1 then Equation (29) in block 210 could be simplified to Q k > n = I k , n and R k . n = 0. 
Additionally if J k - = 1, then accumulator A k . need not be defined and Equation (26), Equation 
(31), and blocks 230 and 240 may be skipped for k=k\ and the summation over k in Equation 
(40) does not include k=k\ As another example, if Y=l (i.e., if Z n = X n in Equation (23)), then 
Equation (34) in block 250 could be simplified to Q n = X n and R n = 0. Additionally if Y = 1 , 
then accumulator B need not be defined and Equation (27), Equation (36), and blocks 270 and 
280 may be skipped, and division by Y in Equation (40) may be skipped. As another example 
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some or all of the intitializations block 200, or mathematical equivalents thereof, could 
alternatively be performed within the iteration loop of FIG. 7, but only at the beginning of the 
first iteration, since the variables Z, B, and A k (k=l, 2. K) are iteration-dependent variables. 

FIG. 8 illustrates use of the algorithm of FIG. 7 to compute Z for the following example: 
Z n = [(51/10)*3 + (42/4)]/6 
and N=5 (i.e., 5 iterations). In this example, X n = (51/10)*3 + (42/4), Y=6, K=2, I ln =51, J,=10, 
M In =3, I 2n =42, J 2 =4, and M 2n =l. A floating point calculation yields X n =25.8, Z n =4.30, Z=21.50 
(i.e., 5*4.30). Thus Z=21.50 is an exact value. 

FIG. 8 shows that the method of FIG. 7, together with application of Equation (40) for the 
post-processing of block 320 of FIG. 7, yields the exact value of Z=21.50. In applying Equation 
(40) and using the values of A,, A 2 , and B from FIG. 8 at the end of iteration 5: 

Z = Z + [B + (A,/J,) + (A 2 /J 2 )]/Y 

= 21 + [2 + (5/10) + (2/4)]/6 
= 21.50 

which agrees with the exact value computed supra. 

FIG. 9 illustrates a computer system 90, in accordance with embodiments of the present 
invention. The computer system 90 comprises a processor 91, an input device 92 coupled to the 
processor 91, an output device 93 coupled to the processor 91, and memory devices 94 and 95 
each coupled to the processor 91. The input device 92 may be, inter alia, a keyboard, a mouse, 
etc. The output device 93 may be, inter alia, a printer, a plotter, a computer screen, a magnetic 
tape, a removable hard disk, a floppy disk, etc. The memory devices 94 and 95 may be, inter 
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alia, a hard disk, a dynamic random access memory (DRAM), a read-only memory (ROM), etc. 
The memory device 95 includes a computer code 97. The computer code or codes 97 includes at 
least one of the algorithms associated with FIGS. 1-3, FIG. 2, FIG. 3, FIG. 6, and FIG. 7, as 
described supra herein. The memory device 94 includes input data 96. The input data 96 
5 includes input required by the computer code or codes 97. The output device 93 displays output 
from the computer code or codes 97. For example, the output device 93 may include 
BA_FINAL output from the algorithm of FIG. 1, BA_FINAL1 output from the algorithm of 
FIG. 2, BAFINAL2 output from the algorithm of FIG. 3, BA output from the algorithm of FIG. 
6, and Z output from the algorithm of FIG. 7. 
10yi While FIG. 9 shows the computer system 90 as a particular configuration of hardware 

y 

jS and software, any configuration of hardware and software, as would be known to a person of 
s ordinary skill in the art, may be utilized for the purposes stated supra in conjunction with the 

o 

In particular computer system 90 of FIG. 9. For example, the memory devices 94 and 95 may be 
portions of a single memory device rather than separate memory devices. 
13* While embodiments of the present invention have been described herein for purposes of 

illustration, many modifications and changes will become apparent to those skilled in the art. 
Accordingly, the appended claims are intended to encompass all such modifications and changes 
as fall within the true spirit and scope of this invention. 
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